Spatial line nodes and fractional vortex pairs in the Fulde-Ferrell-Larkin-Ovchinnikov 

phase 
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A Zeeman magnetic field can induce a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase in spin- 
singlet superconductors. Here we argue that there is a non-trivial solution for the FFLO vortex phase 
that exists near the upper critical field in which the wavefunction has only spatial line nodes that 
form intricate and unusual three-dimensional structures. These structures include a crisscrossing 
lattice of two sets of non-parallel line nodes. We show that these solutions arise from the decay of 
conventional Abrikosov vortices into pairs of fractional vortices. We propose that neutron scattering 
studies can observe these fractional vortex pairs through the observation of a lattice of 1/2 flux 
quanta vortices. We also consider related phases in non-centrosymmetric (NC) superconductors. 

PACS numbers: 



A FFLO phase predicted in Refs. [3,0] appears to have 
been discovered in CeCoIns in the high magnetic field re- 
gion of the superconducting phase diagram |, H- This 
discovery has generated tremendous interest both exper- 
imentally and theoretically ||. FFLO phases have also 
been argued to be of importance in understanding ul- 
tracold atomic Fermi gases Q and in the formation of 
color superconductivity in high density quark matter Q • 
The understanding of these phases has become a rele- 
vant and topical pursuit in physics. One central issue is 
the role vortices play in these phases: in CeCohis the 
FFLO phase appears deep within a vortex phase [!, 0] ; 
and ultracold atomic Fermi gases can be rotated to create 
vortices within an FFLO phase Q. 

Here we address the nature of the FFLO vortex phase. 
Previous studies have concluded that the superconduct- 
ing gap function in this phase is, for example, A(iJ) = 
cos{qz)4> n {r) where the magnetic field is applied along 
the z direction, z ■ r = 0, and <j) n {r) describes a vortex 
lattice constructed from a Landau level (LL) with index 



n [t| [l(J U, 12, HI. This solution has intersecting spa- 
tial nodes along planes perpendicular to the z-axis and 
along the vortex lines parallel to the z-axis. We show 
that there is another realistic solution for the FFLO vor- 
tex phase in which there are only spatial line nodes in 
the gap function. We show that the existence of this 
solution is a consequence of the decay of conventional 
vortices into pairs of fractional vortices. These fractional 
vortices exist because of the broken translational sym- 
metry inherent in FFLO superconductors. By suitably 
choosing an order parameter that correctly exhibits this 
broken translational symmetry, these fractional vortices 
naturally appear within the theory. We propose that 
a small angle neutron scattering (SANS) measurement 
of the resulting magnetic field distribution may observe 
a lattice of 1/2 flux quanta near to the upper critical 
field. We further argue that this phase is stable within 
weak-coupling theories of superconductivity and consider 
related phases in NC superconductors. 

We use a phenomenological approach pioneered by 



Buzdin and Kachkachi to describe the FFLO phase 
U> I2j 13, 14 1, and extend it to include NC supercon- 
ductors. We begin with the following free energy 



F 



J d 3 R{a\A\ 2 + (3\A\ 4 + v\ A| 6 + k|DA| 2 + <5|.D 2 A| 2 
+M|A| 2 |£>A| 2 + 7?[(A*) 2 (£>A) 2 + (A) 2 (£>*A*) 2 ] 

+eB- [A*(DA) + A(DA)*]} (1) 



where D = — iV — 2eA and B = V x A. The coefficients 
that appear in this free energy are typically determined 
from a microscopic BCS theory [3]. The e-term applies 
only to NC superconductors. It results in the helical 
phase discussed previously [lij]. In this phase, the gap 
function becomes A (J?) = ipie^'^. The orientation of 
q is determined by the free energy invariant denoted by 
e in Eq. Q] We have chosen this invariant so that the 
theory applies to Li2Pt3B with point group O [16]. Con- 
sequently, q is parallel to B. With e = 0, Eq. [T]has been 
justified previously [3]. 

We consider a magnetic field along the z direction and 
ignore screening currents in determining the high field 
ground state structure of the gap function (this is rea- 
sonable for strongly type II superconductors). In the 
normal state there will be translational invariance along 
the magnetic field direction. Therefore Fourier modes 
along this direction will be eigenstates of the linear gap 
equation. Typically, the eigenstate with the lowest en- 
ergy corresponds to the Fourier mode q = 0. However in 
FFLO superconductors, the eigenstate with the lowest 
energy have finite q. The states ±g are degenerate and 
this degeneracy is broken by non-linear terms in the free 
energy. Consequently, to describe the FFLO phase near 
the upper critical field, it suffices to keep the two modes 
±q. We therefore write A(R) = vb x (r)e l i z + ip 2 (r)e~ lc t z 
where r is orthogonal to the magnetic field and q is par- 
allel to the field. This yields the following free energy for 
the new order parameter ij> = (ipi, ip 2 ): 

F = L z J dMc^il 2 + a 2 |V> 2 | 2 + /?iM 4 + fo\1>l\ 2 \H 2 

+ 6H</f l^il 2 !^! 2 
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+ Kl (|r>^i| 2 + ir^l 2 ) + K2{\D 2 M 2 + \d 2 H 2 ) 

+ / i(| J D^ 1 | 2 |^i| 2 + |^2| 2 |^| 2 ) 
+ V [(D^(r i ) 2 + {D^ 2 f{r 2 ) 2 + cc] 
+4 V [(D^ 1 )^1(D^ 2 )^ + cc] 
+A4( J DVi)</4(-DV'2)*^2 + c.c] 

rt\DM 2 m 2 + m2\ 2 \M 2 ]} (2) 

where L z is the size of the system along the z direction, 
D = (D x ,D y ) and cc. means complex conjugate. The 
coefficients in Eq. [2] now depend upon q ljj . For FFLO 
superconductors ot\ = a 2 . Eq. [2] should be optimized 
with respect to q and we assume this has been done. 
This ensures that there is no net current flowing along 
the z direction 



ll|,ll5|. 



The choice of order parameter ifi manifestly exhibits 
the broken translational symmetry that characterizes the 
FFLO state. This broken symmetry is hidden when con- 
sidering A. By considering if} explicitly, new and gen- 
eral features of the theory appear naturally. In partic- 
ular, notice that Eq. [2] is independent of separate ro- 
tations of the phases of tpi and tp 2 , revealing a global 
U(l) x U(l) gauge invariance. This follows from transla- 
tional invariance of the normal state along the z direction 
and usual gauge invariance. In particular, consider a gen- 
eral term tp"tp 2 rl (tpl) p (tfj 2 ) q appearing in the free energy, 
usual gauge invariance requires n + m — p — q = and 
translational invariance requires n — m—p+q= 0. These 
two conditions imply that n = p and m = q which leads 
to the U{1) x U{\) invariance. A U(l) x U(l) symmetry 
has been examined to discuss possible topological struc- 
tures in two-band superconductors [lU • Related topolog- 
ical structures have also been discussed in other contexts 

USE!. 

The vortices of a U(l) x U(l) theory can be classified 
(l8| by two integers {n, m) which denote a 2mr phase 
change in i[>i(r) and a 2mir phase change in ip 2 (r) as the 
vortex core is encircled. Of particular interest here are 
the (1,1), (1,0), and (0,1) vortices. The (1,1) vortex is 
the usual Abrikosov vortex and it contains a magnetic 
flux of $o (the usual flux quantum). In the FFLO phase, 
when | = 1 (often called the LO phase), the corre- 
sponding (1, 0) vortex contains a fractional flux $o/2 [H[. 
We are interested in the appearance of bound pairs of 
these vortices in the vortex lattice phase. Consequently, 
we consider generalized Abrikosov vortex lattice states 
and show the usual FFLO vortex solution is often unsta- 
ble to a new lattice solution. In this new solution each 
of the conventional (1,1) vortices decays into a pair of 
(1,0) and (0, 1) vortices. 

We now turn to an analysis valid near the upper crit- 
ical field. The vortex solutions are eigenstates of the 
operator D 2 = (— iV — 2eA) 2 which has eigenvalues 
(2n + l)// 2 and I 2 = <P Q /(2irH) and n = 0,1,2,.. is 
the LL index. The usual BCS theory predicts a n = 
LL solution is the most stable solution, but it has been 



shown that for FFLO superconductors n > LL solu- 
tions can also be stable [5(. It is well known that the LL 
exhibit a macroscopic degeneracy. Abrikosov exploited 
this degeneracy to construct a vortex lattice solution 



which we label as n (^) [221 ] . We label the unit cell of 



the vortex lattice by the lattice vectors a = (a, 0) and 
b = (6 cos a, b sin a). We take r, a, and b to be in units 
Then ab sin a — 2ix gives one flux quantum per unit cell. 
In this basis, we set ip{r) — [r]i(j) n (r) , rj 2 (j) n {r + t)] where 
4> n (r + t) = e~ lTyX '4> n (r + t). The additional phase fac- 
tor that appears in 4> n ensures that both ipi and ip 2 lie 
in the same LL. It appears as a consequence of applying 
a translation in a uniform magnetic field. The new fea- 
ture in this analysis is the appearance of the translation 
vector t = {t Xi t v ) that displaces the nodes of the two 
components , j) 2 ) . Previous results can be recovered 
with r = |lll. |l2j, ll3( . A similar solution has been used 
for UPt3 [23j]. Substituting the above solution for tp(r) 
yields the free energy density (here we have considered 
only the n = LL) 



«i|?7i| 2 +c? 2 |ry 2 | 2 +/3i/3A(0)|r 7 | 4 + 
[(2/3! +/? 2 )/3a(t) - 2/?i/? A (0)]|7 7l | 2 h 2 | 2 + v lA { 



-z/[97a(t) - 37a 



r/| 2 hi| 2 M 2 



(3) 



where the coefficients di,d 2 ,/3i, and /3 2 do not depend 
upon the vortex lattice structure 2J] . The vortex lattice 
structure appears entirely in the generalized Abrikosov 
coefficients Pa{t) = 2n f d 2 r|0 o (»")| 2 |</>o(r + r)| 2 and 
7 a(t) = (2tt) 2 J d 2 r\Mr)\ 4 \Mr + t)\ 2 . Using the 
approach of Ref. [1 21 ] yields 



Mr) 



G 



-GV2giGT 



(4) 



and 



1a(t) = e i| -( GxG ')/V G - T e-( G2 + G ' 2+G - G '>/ 2 (5) 

G,G' 

where G = mgi + ng 2 (n,m are any integer), 
gi = ^2-kgx - v /2irp 2 /ay, g 2 = y/2ir/ay, and 
p + ia = e la b/a. Below, the ground state lattice 
structures are numerically found by minimizing Eq. [3] 
with respect p,cr, and r. 

Single-q to multiple-q transition in NC superconductors. 
Here, c?i ^ d 2 . When di < and d 2 > 0, ryi ^ and 
r/ 2 = 0, the stable structure is the usual hexagonal vortex 
lattice. If/3 2 < 2/3i(/3^ O) (0)-/^(t))//^(t) thenasecond 
transition can occur into a state in which both 771 and 
r/2 are non-zero. This transition has been found within 
weak-coupling theories of NC superconductors [2^, |26|, 
271 ]. This phase has two possible solutions. The first 
has t = and remains a conventional hexagonal lattice. 
This occurs when 2(3i + f3 2 < 0. The second solution 
has t = (a + b)/3 and occurs for 2/3i + (3 2 > 0. To 
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FIG. 1: Helical spatial line nodes in the gap for the multiple-g 
phase of NC superconductors. The centers of the helices form 
a 2D hexagonal lattice perpendicular to the field. 



address which of these possibilities occur, we note that 
weak-coupling microscopic studies show that the phase 
diagram contains a line along which j3 2 = 2(| 27 1. This 
implies that the finite r = (a + b)/3 phase is the ground 
state. 

The spatial nodes of A (it) = e lq2 m<t>o{r) + 
e~ lqz r]2(j)o{r + r) are given by \r]x<po(r)\ = \r] 2 <j)o{r + T )\ 
and cos[qz+(9i — 2 )/2] = where 9\ = 0\(r) is the phase 
of 4>o(r) and 9%(r) is the phase of cj)o (r+r). For small r/ 2 , 
these zeroes lie on small circles surrounding each of the 
zeroes of ipi. Around these circles, the phase 0\{r) = <f) 
since we are encircling a vortex core of ipi (here, </> is the 
polar angle of the circle) and 6 2 ~ cnst since we are far 
away from the zeroes of ^{r). Consequently, the nodes 
of A (it) are given by qz = <p/2 + nn + c where c is a 
constant and n is any integer. This describes the equa- 
tion of a helix spiralling about the z direction. This is 
depicted in Fig. 1. As 772 grows, the pitch of the helix 
grows larger. It is possible for two adjacent helices to 
merge for large enough n 2 . This results in a crisscross- 
ing lattice of line nodes like that discussed below in the 
context of the FFLO case. This analysis reveals that 
the (n,m) = (1,1) Abrikosov vortices have each sepa- 
rated into a pair of (1,0) and (0,1) vortices. The (1,0) 
vortices appear where V'i( r ) = and the (0, 1) vortices 
appear where ^(r) = 0. 

Second order transition into the FFLO phase. Here, cti — 
ot.2, v = 0, and there is one second order transition from 
the normal state into the FFLO state. There are three 
possible solutions for this phase. The first has r\ 2 = 
and rji ^ 0, this is the FF (or single-^) state with a 
conventional hexagonal lattice. This phase is stable when 
(2/3i+/3 2 )/3 a (x)-2/3i/3 j4 (0) >0forall-r. One of the other 
two solutions are stable if (2(3 1 +p 2 )PA{r)-2f3 1 f3 A (0) < 
for any r. The second solution corresponds to \r)i\ = |r/ 2 | 



with r = and is the LO (or multiple-g) phase with a 
conventional hexagonal lattice. This state requires 2/3i + 
(3 2 < to be stable. The final state corresponds to |r?i| = 
1 772 1 with a rectangular unit cell for which b/a = \/3 and 
t = (a, 6)/2 when (3 2 = 0. More generally (including 
first order FFLO transitions) we find the same lattice but 
with b/a 7^ y/3. These solutions are stable for 2/3i + (3 2 > 
0. To understand which of these states may be stable 
within microscopic theories, note that the calculations of 
Ref. [lll | imply that there is a line in the phase diagram 
along which (3 2 = in the weak-coupling theory of a clean 
s-wave superconductor with vortices. Near this line, r = 
(a,b)/2 gives the stable phase. Whenever the FF phase is 
close in energy to the LO phase (that is |/3 2 | « $1), then 
the LO vortex phase with t = (a, b)/2 is the stable vortex 
phase since (3a (t) < /3a (0) for any t^0. It appears that 
this is generic for weak-coupling theories where varying 
gap symmetry, impurities, and vortices lead to a variety 
of different phase diag rams containing both the FF and 
LO phases [n, (H M, El • 

We now focus on the LO phase with r = (a, b)/2. 
This phase can be understood as having conventional 
(1,1) Abrikosov vortices that have each separated into 
a pair of (1,0) and (0,1) vortices. As discussed pre- 
viously the (1,0) and (0,1) vortices in this LO phase 
can be interpreted as containing flux $o/2- To under- 
stand if this may manifest itself experimentally, we have 
performed an Abrikosov analysis [22] on Eq. [2] to deter- 
mine the field distribution h s (r)z due to screening cur- 
rents to lowest order in the gap function. This results in 
h s (r) oc |i/)i(r)| 2 + l^^)! 2 - Consequently, when /3 2 = 0, 
h s (r) has a hexagonal symmetry even though the nodes 
of ipi (r) and ip 2 (r) separately form a rectangular lattice 
with b/a = \^3. A measurement of the hexagonal unit 
cell lattice vector will yield a flux per unit cell that is 
$o/2 (this generalizes to non-hexagonal unit cell geome- 
tries). This can be seen through SANS measurements by 
observing the Bragg peaks of the vortex lattice with neu- 
trons that have momenta perpendicular to the applied 
field. We emphasize that our solution is valid at H c2 and 
at lower fields it is possible that the $o/2 vortices are 
more tightly bound (e.g. r ^ but \t\ < \a + b\/2). 

Here we give the positions of the line nodes for r = 
(a,b)/2 and b/a > 3. In the x,y plane the point zeros 
lie along the lines yi — —36/4, y 2 = —6/4, y 3 = 6/4, 
and 7/4 = 36/4 (the unit cell has doubled along the y- 
direction). The x, z coordinates (measured in units a, ir/q 
respectively) for these four lines are given by z\ = n + 
1/2 — x%/2, z 2 = n + 1/2 + X2/2, Z3 = n — X3/2, and 
Z4 = n + X4/2. This results in a lattice of crisscrossing 
nodal lines as viewed from a direction normal to the y- 
axis (Fig. [2]). 

For the FFLO phase in CeCoIn 5 , the r = (a, 6)/2 so- 
lutions may help to understand some experiments [5(] . In 
particular, measurements in the FFLO vortex phase find 
that the thermal conductivity parallel to the applied field 
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FIG. 2: (Color online) Crisscrossing lattice of nodal lines in 
the FFLO vortex phase with r = (a, b)/2. The dark circles in 
the top figure shows the nodes perpendicular to the applied 
magnetic field for z — 0. As z is changed slightly, these nodes 
move as illustrated by the arrows in the upper right of this 
figure. The lower figure shows a cross section as seen from 
the y direction. The different colors lines correspond to the 
nodal lines with different y positions. 



is g reater than that perpendicular to the applied field 
[29j |. This is not expected for a gap function with spa- 
tial plane nodes perpendicular to the field (which occurs 
if t = 0). However, it can be qualitatively understood 
if t = (a, b)/2. Note that magnetic order in the FFLO 
nodal planes has been proposed (3pj and this may account 
for the thermal conductivity results when r = 0. 

In conclusion, we have argued that the vortex lattice 
phases in FFLO and NC superconductors contain gap 
functions with spatial line nodes that form a variety of 
three dimensional spatial configurations. These configu- 
rations include a lattices of helices in NC superconduc- 
tors and a crisscrossing lattice of nodal lines in FFLO 
superconductors. These structures stem from the break 
up of conventional vortices into pairs of fractional vor- 
tices. SANS studies of the magnetic field distribution 
can provide evidence for these structures. 

We thank Cigdem Capan, Manfred Sigrist, Ilya 
Vekhter, Anton Voronstov, and Kun Yang for useful dis- 



[1] A. I. Larkin and Y. N. Ovchinnikov, Sov. Phys. JETP 
20, 762 (1965). 



[2 

[3; 

[4] 
[5] 

[6] 

[7] 

[8 

[9 

[10 
[11 

[12 

[13 

[i4; 
[is; 

[16 
[17 

[is; 

[19 

[20 

[21 

[22 
[23 
[24 

[25 

[26 

[27 

[28 

[29 
[30 



P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964). 
H.A. Radovan et ai, Nature 425, 51 (2003). 
A. Bianchi et al, Phys. Rev. Lett. 91, 187004 (2003). 
Y. Matsuda and H. Shimahara, J. Phys. Soc. Jpn. 76, 
051005 (2007). 

M.W. Zwierlein et al, Science 311, 492 (2006), G.B. Par- 
tridge et ai, ibid. 311, 503(E) (2006). 
R. Casalbuoni and G. Nardulli, Rev. Mod. Phys. 76, 263 
(2004). 

Y.P. Shim, R.A. Duine, and A.H. MacDonald, Phys. Rev. 
A 74, 053602 (2006). 

T. Mizushima, K. Machida, and M. Ichioka, Phys. Rev. 
Lett. 95, 117003 (2005). 

M. Tachiki et al, Z. Phys. B 100, 369 (1996). 

M. Houzet and A. Buzdin, Phys. Rev. B 63, 184521 

(2001). 

K. Yang and A.H. MacDonald, Phys. Rev. B 70, 094512 
(2004). 

M. Houzet and V.P. Mineev, Phys. Rev. B 74, 144522 

(2006) . 

A. I. Buzdin and H. Kachkachi, Phys. Lett. A 225, 341 
(1997). 

R.P. Kaur, D.F. Agterberg, and M. Sigrist, Phys. Rev. 
Lett. 94, 137002 (2005). 

H. Yuan et al, Phys. Rev. Lett. 97, 017006 (2006). 

Specifically, ot\ — a + nq 2 + 5q 4 + 2eqB, cti — a + nq 2 + 

5q 4 - 2eqB, Ki = k + 26 q 2 , k 2 = S, ft = (3 + fiq 2 + 2riq 2 , 

and ft = 2/3 - 2^q 2 - Ur/q 2 . 

E. Babaev, Phys. Rev. Lett. 89, 067001 (2002). 

M. Sigrist, T. M. Rice, and K. Ueda, Phys. Rev. Lett. 

63, 1727 (1989). 

M.M. Salomaa and G.E. Volovik, Phys. Rev. Lett. 55, 
1184 (1985). 

O. Dimitrova and M.V. Feigel'man, Phys. Rev. B 76, 
014522 (2007). 

A. A. Abrikosov, Soviet Physics JEPT 5, 1174 (1959). 

A. Garg and D.C. Chen, Phys. Rev. B 49, 479 (1994). 
Specifically, &i = ai + ki/I 2 + K2/I 4 , Q2 = a?2 + k,i/1 2 + 
k 2 /1 4 , ft = ft + ii/2l 2 , and ft = ft + ii/l 2 - 

V. Barzykin and L.P. Gor'kov, Phys. Rev. Lett. 89, 
227002 (2002). 

O. V. Dimitrova and M. V. Feigelman, Soviet Physics 
JETP Letters 78, 637 (2003). 

D.F. Agterberg and R.P. Kaur, Phys. Rev. B 75, 064511 

(2007) . 

D.F. Agterberg and K. Yang, J. Phys.: Condens. Matter 
13, 9259 (2001). 

C. Capan et al, Phys. Rev. B 70, 134513 (2004). 

B. L. Young et al, Phys. Rev. Lett. 98, 036402 (2007). 



